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WASP-8b: CHARACTERIZATION OF A COOL AND ECCENTRIC EXOPLANET WITH SPITZER 
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ABSTRACT 

WASP-8b has 2.18 times Jupiter's mass and is on an eccentric (e = 0.31) 8.16-day orbit. With a time-averaged 
equilibrium temperature of 948 K, it is one of the least-irradiated hot Jupiters observed with the Spitzer Space 
Telescope. We have analyzed six photometric light curves of WASP-8b during secondary eclipse observed 
in the 3.6, 4.5, and 8.0 i^m Infrared Array Camera bands. The eclipse depths are 0.113±0.018%, 0.069± 
0.007%, and 0.093 ±0.023%, respectively, giving respective brightness temperatures of 1552, 1131, and 938 
K. We characterized the atmospheric thermal profile and composition of the planet using a line-by-line radiative 
transfer code and a Markov-chain Monte Carlo sampler The data indicated no thermal inversion, independently 
of any assumption about chemical composition. We noted an anomalously high 3.6-^m brightness temperature 
(1552 K); by modeling the eccentricity-caused thermal variation, we found that this temperature is plausible 
for radiative time scales less than ^ 10^ hours. However, as no model spectra fit all three data points well, the 
temperature discrepancy remains as an open question. 

Subject headings: planetary systems — stars: individual: WASP-8 — techniques: photometric 



1. INTRODUCTION 

When transiting exoplanets pass behind their host stars (a 
secondary eclipse), the observed flux drop provides a direct 
measurement of the planet's thermal emission and reflected 
light. Today, secondary-eclipse observations exist for nearly 
30 exoplanets. The Spitzer Space Telescope (Werner et al. 
2004) made most of these observations, capturing broadband 
photometric light curves in six near- and mid-infrared bands 
(3-24 i^m). Each band probes a specific altitude range in a 
planet's atmosphere. With Bayesian fitting of model spec- 
tra, one can quantitatively constrain the atmospheric chem- 
ical composition and thermal profile of the planet's photo- 
sphere (Madhusudhan & Seager 2010). WASP-8b, with a 
time-averaged equilibrium temperature of 948 ± 22 K {T^^, 
temperature at which blackbody emission balances absorbed 
energy, assuming zero albedo and efficient heat redistribu- 
tion), is one of the coolest Jupiter-sized planets yet observed 
in eclipse, and thus serves as an end member to the set of 
measured hot- Jupiter atmospheres. 

To classify the hot- Jupiter population, Fortney et al. (2008) 
proposed a separation between moderately and strongly irra- 
diated planets. The higher atmospheric temperatures of the 
more strongly irradiated planets allow the presence of highly 
opaque molecules (like TiO and VO) at high altitudes. These 
strong absorbers produce hot stratospheres (thermal inversion 
layers). In contrast, for the moderately irradiated hot Jupiters, 
these absorbers condense and rain out to altitudes below the 
photosphere, thus presenting no thermal inversions. 

In general, the observations agree with this hypothesis, but 
exceptions indicate that the picture is not yet completely un- 
derstood. For example, secondary-eclipse observations of 
flie highly irradiated WASP- 12b (Madhusudhan et al. 2011; 
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Crossfield et al. 2012), WASP-14b (Blecic et al. 2013), and 
TrES-3 (Fressin et al. 2010) do not show evidence of ther- 
mal inversions. Conversely, XO-1 has an inversion layer even 
though it receives a much lower stellar irradiation (Machalek 
et al. 2008). Photochemistry provides one explanation. The 
non-equilibrium atmospheric chemistry models of Zahnle 
et al. (2009) suggested that heating from sulfur compounds 
in the upper atmospheres of hot Jupiters could explain these 
inversions. Alternatively, Knutson et al. (2010) suggest that 
strong UV radiation from active stars destroys the high- 
altitude absorbers. 

The Wide-Angle Search for Planets (WASP) Consortium 
discovered WASP-8b in 2008 (Queloz et al. 2010). The planet 
orbits the brighter component (WASP-8A) of a binary stellar 
system. The angular separation (4.83") with the secondary 
(WASP-8B) sets a minimum separation of 440 AU between 
the stars. WASP-8A is a G6 star, with effective temperature 
^eff = 5600 K. Color and photometric analyses indicate that 
WASP-8B is a colder M star (Queloz et al. 2010). WASP-8b 
is a 2.18 Jupiter-mass (Mj^p) planet with 1.08 times Jupiter's 
radius (^j^p) in a retrograde 8.16 day orbit. Its large eccentric- 
ity (e = 0.31) should make the planet's dayside temperature 
vary by hundreds of degrees along the orbit, possibly forcing 
an unusual climate. 

The age of the host star (4 Gyr) is shorter than the planet's 
orbital circularization time (~30 Gyr, see, e.g., Goldreich & 
Soter 1966; Bodenheimer et al. 2001); accordingly, WASP- 
8b has one of the most eccentric orbits among the ~10-day- 
period exoplanets (Pont et al. 2011). The Kozai mechanism 
(Wu & Murray 2003) may explain the combination of high ec- 
centricity and retrograde orbit orientation. The radial- velocity 
(RV) drift and the large eccentricity may also indicate a sec- 
ond planetary companion (Queloz et al. 2010). 

We obtained six secondary-eclipse light curves of WASP- 
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Fig. 1 . — Spitzer images of WASP-8 at 3.6, 4.5, and 8.0 \im, respectively. 
The brighter star, WASP-8A, is at the oiigin. The dimmer WASP-8B signal 
overlapped that of WASP-8A. 

8b from four visits of the Spitzer Space Telescope, observ- 
ing in the 3.6, 4.5, and 8.0 |xm bands of the Infrared Array 
Camera (IRAC, Fazio et al. 2004). The ecUpse depths deter- 
mine the planet's dayside infrared emission. Our Markov- 
chain Monte Carlo (MCMC)-driven radiative-transfer code 
constrained the molecular abundances and temperature pro- 
file of WASP-8b's dayside atmosphere, testing for the ex- 
pected absence of a thermal inversion and estimating the en- 
ergy redistribution over its surface. We constrained the orbit 
of WASP-8b by determining the eclipse epochs and durations. 
We also modeled the thermal variations along the orbit of the 
planet to explore the effects of the eccentricity. 

Section 2 presents the Spitzer observations of the WASP- 
8 system. Section 3 describes the photometric and model- 
ing analysis of our secondary eclipse observations. Section 4 
gives the orbital dynamical analysis. Section 5 presents our 
constraints on WASP-8b's atmospheric composition derived 
from the photometry. Section 6 discusses the effects of eccen- 
tricity on the orbital thermal variation of WASP-8b. Finally, 
Section 7 states our conclusions. 

2. OBSERVATIONS 

The Spitzer Space Telescope visited WASP-8 four times. 
From two consecutive eclipse observations, we obtained si- 
multaneous light curves at 4.5 and 8.0 ixm. Later, from 
two more consecutive eclipse observations during the Warm 
Spitzer mission, we obtained one light curve at 3.6 |J.m and 
one at 4.5 i_im (see Table 1). The Spitzer pipeline (version 
18.18.0) processed the raw data, producing Basic Calibrated 
Data (BCD). 



TABLE 1 
Observation Information 



Label" 


Wavel. 


Observation 


Duration 


Exp. time 


Cadence 




|j.m 


date 


minutes 


seconds 


seconds 


wa()08bs22 


4.5 


2008-12-13 


226 


1.20 


2.0 


wa008bs42 


8.0 


2008-12-13 


226 


10.40 


12.0 


wa008bs21 


4.5 


2008-12-21 


226 


1.20 


2.0 


wa008bs41 


8.0 


2008-12-21 


226 


10.40 


12.0 


waOOSbsU 


3.6 


2010-07-23 


458 


0.36 


0.4 


wa008bs23 


4.5 


2010-07-31 


458 


0.36 


0.4 



" waOOSb designates the planet, i specifies secondary eclipse, and the two 
numbers indicate the wavelength channel and observation serial number (we 
analyzed the 2008-12-21 data before the 2008-12-13 data and inadvertently 
inverted the serial numbers). 

During the initial minutes of our observations, the telescope 
pointing drifted ^--^0.25 pixels before stabilizing. Through- 
out the observations, the pointing also jittered from frame to 
frame (^ 0.01 pixel) and oscillated in an hour-long periodic 
movement (-^ 0.1 pixel amplitude). 



The separation between the centers of WASP-8A and 
WASP-8B in the IRAC detectors is only 3.7 pixels. Con- 
sequently, the signal from the stars overlapped, demanding 
special care during the data analysis (see Figure 1). Table 
2 shows the average and standard deviation of the flux ratio, 
separation, and position angle (PA) of the secondary star with 
respect to WASP-8A (derived from our centering routine, see 
Section 3.2). Our PA values agree with those of Queloz et al. 
(2010), but our separation values are consistently lower than 
theirs (4.83" ±0.01"). 

3. DATA ANALYSIS 

Our Photometry for Orbits, Eclipses, and Transits (POET) 
pipeline produces light curves from BCD images. Briefly, 
POET creates a bad pixel mask for each image, finds the cen- 
ter position of the target, executes interpolated aperture pho- 
tometry, and fits a light curve model that includes physical 
and systematic parameters. 

3.1. POET: Initial Reduction 

POET created bad pixel masks by discarding the flagged 
pixels from the Spitzer BCD masks. Then, it discarded outlier 
pixels with a sigma-rejection method. At each pixel position 
and in sets of 64 consecutive images, POET calculated the 
median and standard deviation of the unmasked pixels. Pixels 
diverging more than four times the standard deviation from 
the median were masked. We iterated this process twice. 

We obtained the Julian Date of each frame from the 
UTCS_OBS and FRAMTIME entries of the files' headers. 
We calculated the Barycentric Julian Date (BJD) by correct- 
ing the projected light- travel time from the telescope to the 
Solar System's barycenter using the Jet Propulsion Labora- 
tory (JPL) Horizons system. We report the times in both 
Coordinated Universal Time (UTC) and the Barycentric Dy- 
namical Time (TDB); the latter is unaffected by leap seconds 
(Eastman et al. 2010). 

3.2. Centering 

POET provides three routines to determine the center of the 
point-spread function (PSF) in each image: center of light, 
2D-Gaussian fitting, and least asymmetry (Stevenson et al. 
2010, Supplementary Information). The proximity of WASP- 
8B confuses these methods, so we added a double-PSF fit 
that shifts supersampled PSFs to the target and secondary, 
bins them down, and scales their amplitudes, as in Crossfield 
et al. (2010). For each Spitzer band we used Tiny Tim' (ver- 
sion 2.0) to create a stellar PSF model with a 5600 K black- 
body spectrum at 100 x finer resolution than our images. The 
double-PSF routine has seven free parameters: the position of 
each star (x^, ja, xb, ys), the integrated stellar fluxes (Fa, Fb), 
and the background sky flux (/sky)- 

To avoid interpolation when binning down, the PSF shifts 
are quantized at the model's resolution, such that image and 
model pixel boundaries coincide. This quantization sets the 
position precision to 0.01 pixels. It also excludes the posi- 
tion parameters from minimizers that assume a continuous 
function, such as Levenberg-Marquardt. So, we fit Fa, Fb, and 
/,ky for a given position set x = {xa, Ja, xb, y^}. 

To avoid the computational challenge of performing a 
minimization for each jc in a 4D space at 0.01 pixel resolution, 

' http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanaly- 
sistools/tools/contributed/general/stinytim/ 
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TABLE 2 
WASP-8 System 



Event 


Flux ratio 


Separation (pixels) 


Separation (") 


Position angle (deg) 




average 


stddev 


average 


stddev 


average 


stddev 


average 


stddev 


waOOSbsU 


0.1420 


0.0030 


3.760 


0.013 


4.610 


0.016 


171.32 


0.28 


wa008bs21 


0.1512 


0.0017 


3.734 


0.007 


4.541 


0.009 


171.00 


0.15 


wa008bs22 


0.1600 


0.0022 


3.737 


0.009 


4.544 


0.011 


170.75 


0.16 


wa008bs23 


0.1648 


0.0039 


3.726 


0.017 


4.497 


0.020 


170.78 


0.29 


wa008bs41 


0.1718 


0.0020 


3.690 


0.009 


4.513 


0.011 


170.84 


0.17 


wa008bs42 


0.1794 


0.0023 


3.686 


0.010 


4.506 


0.012 


170.96 


0.16 



we explored only specific coordinate positions. Starting at 
an initial guess position, and with an initial jump step of 100 
positions ( 1 image pixel), we calculated x' that position and 
the 80 (= 3''- 1) adjacent positions that are one jump step away 
along all combinations of coordinate directions. We either 
moved to the lowest x^or, if already there, shrank the step by 
half. We repeated the procedure until the step was zero. 

3.3. Photometry 

Circular aperture photometry is unsuitable for this system, 
since any flux from the secondary star (WASP-SB) contained 
in the aperture dilutes the eclipse depth of WASP-8b. Small 
pointing jitter would also increase the light-curve dispersion 
for any aperture that included much WASP-8B flux. Aper- 
tures that are too large or small both produce noisier light 
curves. Thus, we modified the POET interpolated aperture 
photometry (Harrington et al. 2007, Supplementary Informa- 
tion) to remove the secondary star two different ways. In both 
methods, we subtracted the median sky level prior to the stel- 
lar flux calculation. The sky annulus included values 7-15 
pixels from the target. 

In our first method (B-Subtract), we subtracted the fitted, 
binned PSF model of WASP-8B from each image. Then, we 
performed interpolated aperture photometry centered on the 
target (A aperture). In the second method (B-Mask), we dis- 
carded the pixels within a circular aperture centered at the po- 
sition of the secondary before performing aperture photom- 
etry. The mask's aperture must encompass most of the con- 
tribution from WASP-8B, but not from the target. Therefore, 
we tested mask apertures with 1.6, 1.8, and 2.0 pixel radii. 
For each photometry method we tested a broad range of A- 
aperture radii in 0.25 pixel intervals. 

The B-Mask method has less residual dispersion when the 
mask is located at a fixed vector separation from WASP- 
8A (using the median of all the measured separations in an 
event), than when its position is determined for each individ- 
ual frame. This can be explained by the dimmer signal of 
WASP-8B, which lowers the accuracy of its position estima- 
tion. So, within each data set using B-Mask, we used the me- 
dian vector separation of the two objects. For the B-subtract 
method, the standard deviation of the normalized residuals 
(SDNR) and eclipse-depths differences are marginal. 

3.4. Light Curve Modeling 

The eclipse depths of WASP-8b are on the order of 0. 1 % of 
the system's flux, well below Spitzer's photometric stability 
criteria (Fazio et al. 2004). Thus, the eclipse light-curve mod- 
eling requires a thorough characterization of the detector sys- 
tematics. Systematic effects have been largely observed and 
documented; they can have both temporal and spatial compo- 
nents, and vary in strength and behavior for each data set. 



The main systematic at 3.6 and 4.5 |J.m is intrapixel sensi- 
tivity variation, M{x,y), where the measured flux depends on 
the precise position of the target on the array (Stevenson et al. 
2012; Charbonneau et al. 2005). In addition, the detectors 
show a time-dependent sensitivity variation called the ramp 
effect, R{t), suspected to be caused by charge trapping (Agol 
et al. 2010) at 8.0 ixm, but there are also reports of a ramp in 
the 3.6 and 4.5 [a.m bands (e.g., Campo et al. 201 1; Nymeyer 
et al. 2011; Knutson et al. 2011; Deming et al. 2011; Blecic 
et al. 2013; Stevenson et al. 2010; Stevenson et al. 2012). The 
eclipse and both systematic variations entangle to produce the 
observed light curve. To account for their contributions, we 
modeled the light curves as 

F(x,y,t) = F,M(x,y)R(t)E(t), (1) 

where Fg is the out-of-eclipse system flux. We used the 
eclipse model, E(t), from Mandel & Agol (2002). The 
eclipse is parametrized by the eclipse depth, the mid-point 
phase, the duration, and the ingress and egress times. For the 
ingress/egress times we adopted a value of 18.8 min, derived 
from the orbital parameters of the planet. We used this value 
in all of our eclipse-model fits. 

The strength and behavior of the ramp variations are spe- 
cific to each dataset. Many formulae have been applied in the 
literature (e.g., Deming et al. 2007; Harrington et al. 2007; 
Knutson et al. 2011; Stevenson et al. 2012). The models are 
formed with combinations of exponential, logarithmic, and 
polynomial functions. We tested dozens of equations; the best 
were: 



risingexp : 


R(t) = 


I - g-n>{t-io) 


(2) 


logramp,^ : 


R(t) = 


\ + r,[\n(t-to)f 


(3) 


linramp : 


R(t) = 


\ + rx{t-tc) 


(4) 


quadrampj : 


R(t) = 


\+rx{t-tc) + r2{t-tcf 


(5) 


quadramp2 : 


R(t) = 


\+r2{t-t,)^ 


(6) 


loglinear : 


R{t) = 


l+ri(f-ff) + r4ln(f-fo) 


(7) 



where tc is a constant value at the approximated mid-point 
phase of the eclipse (tc = 0.5 15 for this planet). Slight changes 
in tc do not significantly affect the fitted eclipse parameters. 

We used our Bi-Linearly Interpolated Subpixel Sensitivity 
(BLISS) mapping technique (Stevenson et al. 2012) to calcu- 
late M{x,y). The BLISS method has been found to return a 
better result than a polynomial fit (Stevenson et al. 2012; Ble- 
cic et al. 2013). 

To determine the best-fitting parameters of our model. 
Equation (1), we used a minimizer with the Levenberg- 
Marquardt algorithm. We used Bayesian posterior sampling 
via an MCMC algorithm to explore the phase space and es- 
timate the uncertainties of the free parameters of the light- 
curve models. Our code implements the MetropoUs random 
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walk, which proposes parameter sets from a muhivariate nor- 
mal distribution centered at the current position in the chain, 
computes x^, and accepts (or rejects) the new set with greater 
probability for a lower (higher) x^. By generating millions 
of parameters sets, the algorithm samples the posterior dis- 
tribution of the model parameters. As a necessary condition 
for chain convergence, we require the Gelman-Rubin statistic 
(Gelman & Rubin 1992) to be within 1% of unity for each 
free parameter between four MCMC chains. 

The photometry routine uses the BCD uncertainty images 
to estimate the uncertainties of the light-curve data points, cr, . 
However, since the Spitzer pipeline in general overestimates 
these uncertainties (it is designed for absolute photometry), 
we multiply by a constant factor (cr, f ■ cr,), such that the 
reduced = 1 in the light-curve fit. This is equivalent to 
estimating a single cr from the scatter of model residuals. Both 
methods account for red noise, but ours retains the (usually 
small) cr, variations due to aberrant frames. 

To determine the best raw light curve (i.e., by selection of 
photometry method and aperture radius), we calculated the 
SDNR of the light-curve fit (Stevenson et al. 2012; Campo 
et al. 2011). Poor fits or data with high dispersion increase 
SDNR; the optimum data set minimizes the SDNR value. 
Once we chose the best light curve, we compared the different 
ramp models according to the Bayesian Information Criterion 
(Liddle 2007), 

BlC = x^ + k\nN, (8) 

where k is the number of free parameters and the num- 
ber of data points. The best model minimizes the BIC. The 
probability ratio favoring one model over a second one is 
exp(- ABIC/2). 

3.4.1. waOOSbsll Analysis 

This observation started 2.9 hours before the eclipse's first 
contact. The telescope observed the target in sub-array mode, 
allowing a high cadence (Table 1). We discarded the initial 15 
minutes of observation while the telescope pointing settled. 
Our data present both intrapixel and weak ramp systematics. 



TABLE 3 
waOOSbs 1 1 Ramp Model Fits 



R(t) 


SDNR 


ABIC 


Eel. Depth (%) 


quadrampi 


0.0061141 


0.00 


0.119 


risingexp 


0.0061148 


1.73 


0.106 


logramp j 


0.0061153 


2.96 


0.096 


linramp 


0.0061201 


6.34 


0.063 


loglinear 


0.0061141 


11.10 


0.119 



The 2.25 pixel A aperture with B-subtract photometry mini- 
mized SDNR. Table 3 shows the five best-fitting models to the 
best waOOSbsl 1 light curve. ABIC is with respect to the low- 
est BIC value. The quadrampj model is 2.4 times more prob- 
able than, and consistent with, the second-best model. The 
linear (llcr) and and quadratic (4cr) terms of the quadramp; 
model (see Table 9) confirm the need for a ramp model. As 
a general remark, we noted that all the logramp ^ models pro- 
duce similar BIC and eclipse parameter values; therefore, we 
will refer only to the logramp j model in the future. Mod- 
els with more free parameters do not improve BIC. Following 
Stevenson et al. (2012), we vary the bin size and the minimum 
number of data points per bin (mnp) of the BLISS map to min- 
imize the dispersion of the residuals. We required at least 4 



points per bin for any dataset. The PSF-fitting position pre- 
cision of 0.01 pixels sets our lower limit for the binsize. For 
waOOSbsl 1, mnp = 5 and a bin size of 0.015 pixels optimized 
the fit. 

Figure 2 shows the raw, binned and systematics-corrected 
waOOSbsll light curves with their best-fitting model. We 
considered the correlated noise in the residuals as well (Pont 
et al. 2006). Figure 3 shows the root-mean-square (RMS) 
of the residuals vs. bin size. The waOOSbsll RMS curve 
deviates above the expected RMS for pure Gaussian noise. 
Following Winn et al. (200S), to account for the correlated 
noise, we weighted the light curve uncertainties by the fac- 
tor f3 (the fractional RMS excess above the pure Gaussian 
RMS at the bin size corresponding to the eclipse duration). 
For waOOSbsl 1 we found /? = 2.4. We inspected all the pair- 
wise correlation plots and histograms and found only uni- 
modal Gaussian distributions. 

Alternatively, the residual-permutation (also known as 
prayer bead) algorithm is sometimes used to assess correlated 
noise in a fit. In this method, we cyclically shift the residu- 
als from the best model by one frame, add them back to that 
model, and re-fit, repeating until we shift the residuals back to 
their original positions. This generates a distribution of values 
for each parameter, from which we estimate the parameter un- 
certainties. The eclipse-depth uncertainty is 0.021%, similar 
to the value found with the Winn et al. (200S) method (see Ta- 
ble 9). However, we are cautious. Although prayer bead has 
been broadly used for the analysis of exoplanet lightcurve and 
RV fits (e.g., Southworth 200S; Bouchy et al. 2005; Pont et al. 
2005; Gillon et al. 2007; Knutson et al. 200S; Cowan et al. 
2012), we have found no detailed description of its statistical 
properties in the literature. 

3.4.2. wa008bs2 3 Analysis 

With the same observing setup as waOOSbsll, this obser- 
vation started 3.3 hours prior to the eclipse's first contact. 
This dataset also presented both intrapixel and ramp vari- 
ations. Note that the intrapixel systematic is weaker than 
at 3.6 \xm, attributed to the smaller degree of undersam- 
pling at larger wavelengths by Charbonneau et al. (2005) and 
Morales-Calderon et al. (2006). Even though the pointing sta- 
bilized only after the initial 20 minutes, the light curve did not 
deviate significantly; therefore, we included all data points in 
the analysis. We noted two sudden pointing and PA deviations 
near phase 0.519. After each incident, the telescope resumed 
its position within ~10 s (Figure 4). Micrometeorite impacts 
on the telescope can explain the abrupt deviations. Simulta- 
neously, we measured a slight increase in the background sky 
flux dispersion, which returned to normality shortly after The 
target flux did not show any extraordinary fluctuations during 
these incidents. However, the points outside the normal point- 
ing range were eliminated by the BLISS map's mnp criterion. 

The analysis is analogous to waOOSbsll. The SDNR indi- 
cated clearly that the 2.25 pixel A aperture with B-Subtract 
photometry produced the lowest dispersion. The best-fitting 
ramp model is logrampj, which is 21 times more probable 
than the rising exponential ramp (Table 4). The BLISS map 
is optimized at mnp = 4 and a bin size of 0.025 pixels. 

An initial MCMC run showed a significant linear correla- 
tion between the system flux and the r^ parameter of the lo- 
gramp, which prevented the MCMC chain from converging. 
We solved this problem by transforming the correlated pa- 
rameters into an orthogonal set of parameters, rerunning the 
MCMC chain, and inverting the transformation on the result- 
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Fig. 2. — Raw (left), binned (center) and systematics-corrected (right) secondary-eclipse light curves of WASP-8b at 3.6, 4.5, and 8.0 |xm. The system flux 
is normalized to unity and the points are shifted vertically for clarity. The colored curves are the best-fit models (see legend). The black curves are the best-fit 
models excluding the eclipse component. The error bars in the center and right panels are the la unceitainties. 
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Fig. 3. — RMS of the fit residuals (black curves with Icr uncertainties) vs. 
bin size of the WASP-8b light curves. The red curves are the expected RMS 
for Gaussian noise (extrapolated unbinned RMS scaling by the inverse square 
root of the bin size). The blue dotted and green dashed vertical lines mark the 
bin size corresponding to the eclipse ingress and duration time, respectively. 
wa008bs 11 's excess above the red line indicate correlated noise at bin sizes 
larger than the ingress time. 

ing parameter values (Stevenson et al. 2012). Figures 2 and 3 
show the wa008bs23 light curves with the best-fitting model 
and RMS of the residuals vs. bin size. 

3.4.3. wa008bs21 cSc wa008bs41 Analysis 

We simultaneously observed wa008bs21 and wa008bs41 
in full-array mode. Prior to the eclipse observation, we ex- 



TABLE4 
wa008bs23 Ramp Model Fits 



R(t) 


SDNR 


ABIC 


Ed. Depth (%) 


logramp j 


0.0073830 


0.00 


0.0677 


risingexp 


0.0073832 


6.10 


0.0730 


quadrampi 


0.0073833 


9.40 


0.0777 


loglinear 


0.0073830 


10.84 


0.0685 


linramp 


0.0073846 


12.19 


0.0564 



posed the detector (a "preflash" observation, Knutson et al. 
2009) for 25 minutes to a bright H II region, with coordinates 
a = 20'' 21'" 39128 and (5 = +37° 31' 03.6", to minimize the 
ramp systematic variation. The secondary-eclipse observation 
started only 26 minutes before the first contact. The telescope 
pointing stabilized quickly, so fortunately we needed to re- 
move only the initial four minutes of observation. Every 12 
s, the detector recorded two consecutive images (two-second 
exposures) at 4.5 |am and one image at 8.0 ^.m. (Table 1). 

The SDNR analysis of wa008bs21 showed that a 3.5 pixel 
A aperture with 1 .6 pixel B-Mask photometry minimizes the 
dispersion (Figure 5). The ramp models indicated a negligi- 
ble ramp variation. Accordingly, a fit without a ramp model 
yielded the lowest BIC. Table 5 shows the four best-fitting 
models for the best wa008bs21 data set. The no-ramp model 
is 15 times more probable than the quadramp2 model. 

Because of the shorter out-of-eclipse observation, the sys- 
tem flux is less-constrained for wa008bs21 than for the 
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Fig. 4. — wa008bs23 target pointing, position angle and background sky 
flux near phase 0.519. We observed two sudden position shifts coincident 
with increases in the background flux. All values returned to normal almost 
instantly. 
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B aperture = 1.6 
B aperture = 1.8 
B aperture = 2.0 
subtract B 




3.00 3.25 3.50 3.75 4.00 4.25 

A aperture 

Fig. 5. — wa008bs21 Standard Deviation of the Normalized Residuals vs. 
aperture (in pixels). The SDNR curves use the best ramp model in Table 5. 
The legend indicates the photometry method. SDNR increases at 3.75 pixels 
(coincident with the stars' separation). The eclipse parameters are consistent 
over the 3.0-4.25 aperture range. The optimum dataset uses 1 .6 pixel B-Mask 
photometry with a 3.5 pixel A aperture. 



TABLE 5 
WA008BS21 Ramp Model Fits 



R(t) 


SDNR 


ABIC 


Eel. Depth (%) 


no ramp 


0.0036223 


0.00 


0.0718 


quadramp2 


0.0036195 


5.76 


0.1189 


linramp 


0.0036223 


7.56 


0.0714 


quadrampi 


0.0036197 


13.14 


0.1170 



waOOSbsll or wa008bs23 events. Combined with a corre- 
lation between the eclipse depth and system flux (revealed 
by MCMC), the lower precision of the system flux trans- 
lates into a larger eclipse depth uncertainty. Nevertheless, the 



wa008bs21 fit parameters were consistent among the different 
apertures (Figure 6). The optimum parameters of the BLISS 
map are mnp = 5 and a bin size of 0.025 pixels. 
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3.00 3.25 3.50 3.75 4.00 4.25 
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Fig. 6. — Eclipse depth V.V. A aperture for wa008bs21. Each color represent 
a different photometry method as in Figure 5. The blue en'or bar corresponds 
to the \-a uncertainty of the best model. The eclipse duration and mid-point 
phase follow a similar trend. 

The 8.0 i^m detector did not present an intrapixel pattern 
like the 3.6 or 4.5 |xm detectors. However, some of the raw 
light curves for different apertures and photometry methods 
showed large scatter and presented strong oscillations, pro- 
ducing implausible fit parameters. A pixelation effect (An- 
derson et al. 2011; Stevenson et al. 2012) might be respon- 
sible. As a consequence, we were unable to fit the eclipse 
parameters unambiguously for this data set alone. Normally 
we study the events individually to select the best aperture and 
photometry method, but in this case we used a joint fit with 
the best wa008bs21 dataset and model to help constrain the 
8.0 ^.m eclipse curve, sharing the eclipse duration and mid- 
point parameters. The 3.5 pixel A aperture with 1.6 pixel B- 
mask photometry for wa008bs41 minimized the joint SDNR 
(Figure 7). 

Table 6 compares the four best-fitting ramp models for the 
best wa008bs41 light curve. A linear ramp minimizes BIC, 
and is 20 times more probable than the next-best model. Fig- 
ures 2 and 3 show the wa008bs21 and wa008bs41 light curves 
with their best-fitting models and RMS of the residuals vs. bin 
size, respectively. 



B aperture = 1.6 
B aperture = 1.8 
B aperture = 2.0 
B Subtract 




3.00 3.25 3.50 3.75 4.00 4.25 

A aperture 

Fig. 7. — Joint wa008bs21+wa008bs41 Standard Deviation of the Normal- 
ized Residuals vs. aperture radius (in pixels) of wa008bs41, for different pho- 
tometry methods. All 24 light curves use the best ramp model from Table 6. 
Light curves using 1.6 and 2.0 pixel B-Mask photometry at 3.5 pixel A aper- 
ture produce consistent eclipse parameters and outpeiform the best B-subtract 
method (with a 3.75 pixel A aperture). The best B-subtract also yields a more 
scattered raw light curve. Hence, the optimum dataset uses 1 .6 pixel B-mask 
photometry with a 3.5 pixel A aperture. 
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TABLE 6 
wa008bs41 Ramp Model Fits 



R(t) 


SDNR 4&2 


ABIC 


Eel. Depth (%) 


linramp 


0.0030766 


0.00 


0.0931 


quadrampi 


0.0030744 


5.94 


0.1308 


risingexp 


0.0030754 


6.33 


0.1150 


logrampi 


0.0030758 


6.50 


0.1078 



3.4.4. wa008bs22 & wa008bs42 Analysis 

The observing setup of these events was identical to 
wa008bs21 and wa008bs41, including the preflash observa- 
tion. The pointing of this observation drifted noticeably more 
than in the other observations, moving more than 0.4 pixels 
during the initial 30 minutes and stabilizing only during the 
eclipse. As a consequence, the illumination level of the indi- 
vidual pixels changed during the beginning of the eclipse. The 
ramp variation, which depends on the illumination (Knutson 
et al. 2008), was disrupted. 

The wa008bs22 event, having a negligible ramp variation, 
was little affected by the telescope pointing shift. The SDNR 
calculation for wa008bs22 indicated the 1.8 pixel B-Mask 
photometry with 3.75 pixel A aperture as the best dataset. A 
light-curve model without a ramp (Table 7) is 639 times more 
probable than the quadramp2 model. The optimal BLISS map 
has mnp = 4 and a bin size of 0.02 pixels. 



TABLE 7 
WA008BS22 Ramp Model Fits 



R(t) 


SDNR 


ABIC 


Eel. Depth (%) 


no ramp 


0.0025274 


0.00 


0.0814 


quadramp2 


0.0025253 


12.94 


0.1224 


logramp j 


0.0025267 


14.22 


0.0921 


risingexp 


0.0025274 


15.04 


0.0814 



In contrast, we discarded the initial wa008bs42 light curve 
past the eclipse ingress due to the disrupted ramp variation. 
The eclipse model parameters are thus less constrained. By 
this point, we already had single-channel fits for the rest of 
the data, so we tuned the wa008bs42 analysis in a joint fit 
with all the other events, sharing the eclipse duration and 
mid-time. SDNR indicates the B-subtract method with 4.00 
pixel A aperture as the best dataset. Table 8 presents the four 
best-fitting models. The eclipse depth is consistent with the 
wa008bs41 depth. 



TABLE 8 
wa008bs42 Ramp Model Fits 



R(t) 


SDNR 


ABIC 


Eel. Depth (%) 


linramp 


0.0032320 


0.00 


0.0932 


quadrampi 


0.0032312 


6.19 


0.0892 


logramp j 


0.0032321 


6.67 


0.0938 


risingexp 


0.0032330 


7.08 


0.0961 



3.4.5. Final Joint-fit Analysis 

From the three individual fits to the 4.5 i_im observations 
we found eclipse depths of 0.072% ± 0.021%, 0.086% ± 
0.022%, and 0.068% ± 0.007% for wa008bs21, wa008bs22, 
and wa008bs23, respectively. The weighted mean of the 
depths is 0.0692% ± 0.0065%. With a dispersion of 0.0062% 



around the mean, this is not larger than the individual un- 
certainties, thus we found no evidence for temporal variabil- 
ity. This dispersion corresponds to 10% of the mean eclipse 
depth. The consistency permitted a joint analysis of all obser- 
vations. We used the best light curves and models found in the 
individual fits, where all events shared the eclipse duration, 
the three 4.5-iam events shared their eclipse depth, the two 
8.0- [a.m events shared their eclipse depth, the simultaneous 
wa008bs21 and wa008bs41 events shared their eclipse mid- 
point phases, and the wa008bs22 and the wa008bs42 events 
shared their eclipse mid-point phases. Table 9 shows the light- 
curve modeling setup and results. We used these joint-fit re- 
sults for the orbital and atmospheric analysis. An electronic 
supplement contains the best light curves, including center- 
ing, photometry, and the joint fit. 

4. ORBITAL DYNAMICS 

WASP-8b's high eccentricity (e = 0.31) implies that its 
separation from WASP-8A at periapsis (0.055 AU) is about 
half that at apoapsis. Given the argument of periapsis (cj = 
-85.86°), the secondary eclipse nearly coincides with the pe- 
riapsis. The planet, therefore, receives over twice as much 
flux at eclipse as it would if the orbit were circular, explaining 
in part our high brightness temperature (see Table 9). 

Secondary-eclipse times can refine estimates of ecoso; 
from RV data. The four eclipse events occurred at an aver- 
age eclipse phase of 0.514695 ± 0.00018. After subtracting 
a coarse light-time correction of 2fl/c = 80 s from this av- 
erage phase, we calculated ecoscj = 0.02290 ± 0.00028 (see 
Equation (3) of Charbonneau et al. 2005). This is consistent 
with Queloz et al. (2010), and photometrically confirms the 
nonzero eccentricity of the planet's orbit (we fit ecosu below 
without relying on the low e approximation). 

The eclipse timings were combined with 130 available RV 
data points and with transit data from Queloz et al. (2010) 
using the method described by Campo et al. (2011) and 
Nymeyer et al. (201 1). Forty-eight in-transit RV points were 
removed due to the Rossiter-McLaughlin effect. 

Our fit presented a moderate improvement to the orbital 
parameters of WASP-8b (Table 10), except for the period. 
While Queloz et al. (2010) used several transits to mea- 
sure the period, we used their published mid-point epoch 
(a single date); hence, our period is constrained mostly by 
our eclipses and the RV data, and thus have a larger un- 
certainty. By themselves, the secondary eclipses have a pe- 
riod of 8.158774 ±0.00040 days and a midpoint epoch of 
BJD 2455409.6629 ± 0.0017 (TDB), not significantly ([5.9 ± 
4.3] X 10"^ days) longer than the period found by Queloz et al. 
(2010). The transit and eclipse periods place a 9.8 x 10"^ 
° day"' (3(t) upper limit on possible apsidal precession, nearly 
three orders of magnitude larger than the theoretical expecta- 
tion for tidal effects (Ragozzine & Wolf 2009). 

5. ATMOSPHERIC ANALYSIS 

We use our IRAC observations of thermal emission from 
WASP-8b to constrain the thermal structure and composition 
of the dayside atmosphere of the planet. The Spitzer band- 
passes at 3.6, 4.5, and 8.0 \xm contain strong spectral features 
due to several carbon and oxygen-based molecules that are ex- 
pected in hot- Jupiter atmospheres. Methane (CH4) has strong 
spectral features in the 3.6 and 8.0 |xm bands, carbon monox- 
ide (CO) and carbon dioxide (CO2) have features at 4.5 |xm, 
while water vapor (H2O) has features in all three bands (Mad- 
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Notes: The values quoted in parenthes are the Icr uncertainties. 
^RMS frame-to-frame position difference. 
''MJD = BJD - 2,450,000. 

■^We exclude frames duiing instrument/telescope settling, for insufficient points at a given BLISS knot, and for bad pixels in the photometry aperture, 
■'in the individual fits. Joint fit had 19 free parameters. 





TABLE 10 






Eccentric Orbital Model 


Parameter 


This work 


Queloz et al. (2010) 


esinci; 


-0.3078 ± 0.0020 


-0.3092 ±0.0029 


ecosw 


0.02219 ± 0.00046 


0.023 ±0.001 


e 


0.309 ±0.002 


0.310 ±0.0029 


aj(°) 


-85.00 ±0.08 


-85.73 ±0.18 


P (days) 


8.158719 ±0.000034 


8.158715 ±0.000016 


To (MJDtdb) 


4679.33486 ± 0.00057 


4679.33509 ± 0.00050 


A:(ms'') 


221.9 ±0.6 


222.23 ±0.8 


7c (ms'') 


-1565.9 ±0.6 


-1565.8 ±0.21 


7h (ms^') 


-1547.4 ±0.4 


-1548.1 ±0.6 


7 (ms"'yr~') 


58.1 ± 1.2 


58.1 ±1.3 


Reduced X" 


4.1 


0.86 



husudhan & Seager 2010). The spectral features of the vari- 
ous molecules appear as absorption troughs or emission peaks 
in the emergent spectrum depending on whether the temper- 
ature decreases or increases with altitude, respectively. Con- 
sequently, strong degeneracies exist between the temperature 
structure and molecular composition derived from a spectral 
dataset (e.g., Madhusudhan & Seager 2010). Nevertheless, 
photometric observations made with Spitzer have been suc- 
cessfully used to constrain chemical compositions and tem- 
perature structures in many exoplanetary atmospheres (e.g.. 
Barman et al. 2005; Burrows et al. 2007; Knutson et al. 2008; 
Madhusudhan & Seager 2009; Stevenson et al. 2010; Mad- 
husudhan et al. 201 1). 



We model the dayside emergent spectrum of WASP-8b us- 
ing the atmospheric modeling and retrieval method of Mad- 
husudhan & Seager (2009, 2C)10). The model computes line- 
by-line radiative transfer in a plane-parallel atmosphere as- 
suming hydrostatic equilibrium, local thermodynamic equi- 
librium, and global energy balance. We assume a Kurucz 
model for the stellar spectrum (Castelli & Kurucz 2004) given 
the stellar parameters. The pressure-temperature (P-T) pro- 
file and molecular mixing ratios are free parameters in the 
model, which can be constrained from the data. The P-T 
profile comprises of six free parameters and the mixing ra- 
tio of each molecular species constitutes an additional free 
parameter. Following Madhusudhan & Seager (2009), we 
parametrize the mixing ratio of each species as deviations 
from thermochemical equilibrium assuming solar elemental 
abundances (Burrows & Sharp 1999). We include the domi- 
nant sources of opacity expected in hot Jupiter atmospheres, 
namely molecular absorption due to H2O, CO, CH4 and CO2 
(Freedman et al. 2008; R.S. Freedman 2009, private commu- 
nication), and H2-H2 collision induced absorption (Borysow 
2002). We explore the model parameter space in a Bayesian 
way using an MCMC sampler (Madhusudhan & Seager 2010, 
2011). Given the limited number of observations (Wobs=3), 
our goal is not to find a unique model fit to the data; instead, 
we intend to constrain the region of atmospheric parameter 
space that is allowed or ruled out by the data. 

Our observations rule out a thermal inversion in the dayside 
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A iiivn) Normalized contribution 



Fig. 8. — Left: Atmospheric spectral emission of the dayside of WASP-8b. The blue circles with error bars are the measured echpse depths, or equivalently, the 
planet-star flux ratios. The red and green curves show two model spectra with the same temperature profile (shown in the inset) but with different compositions. 
The green curve shows a model assuming chemical equilibrium with solar elemental abundances. The red curve shows a model with lO' times lower methane 
abundance compared to the green model, but the abundances of the remaining molecules are identical to those in the green model. The red and green filled circles 
are the corresponding model fluxes integrated over the Spitzer bands (bottom solid lines). The black dashed lines represent planetary blackbody spectra with 
r = 710, 1100, and 1450 K. Right: Normalized contribution functions of the solar-composition model (solid lines) and the low-CH4-abundance model (dotted 
lines) in each Spitzer band (see legend). The effective pressures of the contribution functions are 0.63, 0.35, and 0.12 bar at 3.6, 4.5, and 8.0 i^m, respectively. 



atmosphere of WASP-8b. This is evident from the planet-star 
flux contrasts in the three IRAC bands at 3.6, 4.5, and 8.0 
|j.m. In the presence of a thermal inversion, the planet-star 
flux contrasts in the 4.5 and 8.0 |xm bands are both expected 
to be greater than the flux contrast in the 3.6 |im band (Bur- 
rows et al. 2008; Fortney et al. 2008; Madhusudhan & Seager 
2010), due to spectral features of the dominant molecules ap- 
pearing as emission peaks as opposed to absorption troughs. 
However, the low 4.5 and 8.0 [im flux contrasts relative to the 
3.6 |J.m contrast requires significant absorption due to H2O 
and CO across the spectrum, and hence the lack of a thermal 
inversion in the atmosphere. Figure 8 shows model spectra of 
WASP-8b with no thermal inversion in the temperature pro- 
file. The observed 4.5 and 8.0 ^.m flux contrasts are explained 
to a good level of fit by a model without a thermal inversion 
and with solar abundance composition, as shown by the green 
curve in Figure 8. Our inference of the lack of a thermal in- 
version in WASP-8b is independent of any assumption about 
chemical composition or C/O ratio (e.g. Madhusudhan & Sea- 
ger 201 1). The lack of a thermal inversion in WASP-8b is not 
surprising, since it is amongst the cooler population of irra- 
diated hot Jupiters, which are not expected to host inversion- 
causing species such as TiO or VO in their upper atmosphere 
(Fortney et al. 2008; Spiegel et al. 2009). 

Our models are unable to reproduce the high planet-star 
flux contrast observed in the 3.6 ixm IRAC band, indepen- 
dent of the composition. The major sources of absorption in 
the 3.6 [a.m band are H2O and CH4. In principle, decreasing 
the CH4 and/or H2O abundances can lead to a higher 3.6 i_im 
contrast. However, as shown by the red curve in Figure 8, 
such an increase also simultaneously increases the contrast in 
the 8.0 i_im band, thereby worsening the fit overall. Another 
hindrance to fitting the observed 3.6 [xm contrast is that it re- 
quires a hotter P-T profile, with T > 1550 K in the lower 
atmosphere, predicts much higher fluxes in the 4.5 and 8.0 
|j.m bands than observed. On the other hand, a cooler P-T 



profile than shown in Figure 8 would provide a better fit in 
the 4.5 and 8.0 ^.m bands, but would further worsen the fit in 
the 3.6 i^m band. Consequently, we choose an intermediate 
P-T profile that provides a compromise fit to all three data 
points. 

Although the one-dimensional (ID) models shown in Fig- 
ure 8 output less energy than the instantaneous incident ir- 
radiation during the eclipse (concurrent with periastron pas- 
sage), they output ^ 20% higher energy compared to the time- 
averaged incident irradiation received at the substellar point. 
Considering that a pseudo-synchronous rotation should facil- 
itate the redistribution of energy to the night side, the high 
emission measured suggests that WASP-8b is quickly reradi- 
ating the incident irradiation on its dayside hemisphere, i.e. 
nearly zero day-night redistribution. Such a scenario would 
lead to a large day-night temperature contrast in the planet 
which can be confirmed by thermal phase curves of the planet 
observed using warm Spitzer (e.g., Knutson et al. 2009). The 
high emergent flux also implies a very low albedo, as with 
most hot- Jupiter planets (Cowan & Agol 201 lb). 

6. THE UNEXPECTED BRIGHTNESS TEMPERATURES 
OF WASP-8b 

As seen in the previous section, the 3.6-|^m brightness 
temperature is anomalously higher than expected. The 
hemisphere-averaged equilibrium temperature for instanta- 
neous reradiation (time-averaged around the orbit) is only 948 
K; even the instantaneous equilibrium temperature at periap- 
sis, 1128 K, is far lower than this observation. Thus, we mod- 
eled the orbital thermal variation due to the eccentricity to de- 
termine if such a high temperature is possible from irradiation 
alone. 

Following Cowan & Agol (2011a), we solved the energy 
balance equation in a one-layer latitude-longitude grid over 
the planetary surface. The change in temperature of a cell 
with time, dr/df, is determined by the difference between the 



10 



Cubillos et al. 



absorbed flux from the star and the re-emitted blackbody flux, 

2 



dT 
dT 



1 



R., 



(9) 



where c/, is the heat capacity per unit area; T^ff and are the 
star's effective temperature and radius, respectively; r{t) is the 
planet-star separation; cos = sin Amax(cos$(f),0) is the 
cosine of the angle between the vectors normal to the planet 
surface and the incident radiation, with A the latitude of the 
cell and $(?) the longitude from the substellar meridian; a is 
the Stefan-Boltzmann constant. 

Tidal interactions drive the planet's rotational angular ve- 
locity (Wjot) toward synchronization with the orbital angular 
velocity (worb)- Hence, if the spin synchronization timescale 
(e.g., Seager & Hui 2002; Goldreich & Soter 1966) is shorter 
than the system age, we expect cji-ot = Woit- In the case of 
WASP-8b, the timescale for tidal synchronization is on the 
order of 0.05 Gyr, much shorter than the age of the star How- 
ever, a planet in an eccentric orbit, where Worb changes in 
time, is actually expected to reach a pseudo-synchronization 
state (e.g., Langton & Laughlin 2008; Hut 1981), in which 
the planet does not exchange net angular momentum with its 
orbit. The planet acquires then a constant rotational angu- 
lar velocity close to the orbital angular velocity at periastron 
(worb.p)- In the literature we found different predictions for 
this equilibrium angular velocity, from 0.8a;oi-b,p (Hut 1981) 
to l.SScjoi-b.p (Ivanov & Papaloizou 2007). 

The tidal evolution drives the orbit of a planet toward zero 
obliquity in a timescale similar to the spin synchronization 
(Peale 1999). We thus adopted zero obliquity for our simula- 
tions. We also assumed A = 0, supported by the atmospheric 
analysis (Section 5). Beyond these assumptions, the parame- 
ters of interest that control Equation (9) are the radiative time 
Trad = Ch/(jTQ (whcrc Tq is the substellar equilibrium temper- 
ature at periastron) and the rotational angular velocity of the 
planet Wiot (which determines the substellar longitude of a cell 
through the equation d$(f )/df = uji-ot-LOorbit))- With these def- 
initions Equation (9) can be re-written as: 



dT 
dT 



Trad 



max(cos$(0,0) - ( 
r(t) J \ToJ 



(10) 



We derived the temperature of each cell as a function of 
time to study its thermal evolution. Assuming that each cell 
emits as a blackbody, we calculated the photometric phase 
curve of the planet by integrating over the hemisphere observ- 
able from Earth, weighted by the viewing geometry. Our sim- 
ulations were for planets nearly in pseudo-synchronous rota- 
tion (cjrot = 0.8, 1.0, and 1.5 Worb.p)- We tested values of t,.^^ 
between 1 and lO"* hours. 

Figure 9 shows simulated brightness-temperature 
lightcurves of WASP-8b after reaching a periodic sta- 
tionary state (after a few t^.^^). We noted that the higher 
irradiation at periastron is not the only contribution to a 
higher temperature. For Wrot > Worb,p, the substellar angular 
velocity (d<l> / df ) is minimum during periastron, allowing the 
temperature to increase due to the longer exposure to the 
irradiation. For Wrot < Worb.p, the substellar angular velocity 
is negative for an instant around periastron. Later, when the 
planet emerges from secondary eclipse the over-heated region 
becomes observable from Earth. As a result, the lightcurve 
shows a delayed maximum. 



E 1000 




0.25 0.5 0.75 1.0 

Orbital pbase (8.159-day period) 

Fig. 9. — Model brightness-temperature lightcurves of WASP-8b as ob- 
served from Earth during one orbital period. Phase zero indicates the mid- 
transit time. The gray region indicates the secondary-eclipse interval, with 
periastron at phase 0.52. The models simulate super-rotating winds (Wrot = 
1.5 LOoib.p) for different radiative times (see legend). The curves with smaller 
Trad show larger amplitudes. For t,-^^ comparable to the orbit period, and 
since a;,-ot =1.5 uja^h p, opposite sides will face the star during successive pe- 
riastron passages, leading to two bright spots and hence three periodic peaks 
per orbit. 



Our models show that for large radiative timescales, the 
temperatures at secondary eclipse are lower than 1 150 K, re- 
gardless of ojrot- For radiative times shorter than ^ 10^ hr, the 
temperatures can be as high as 1400 K, similar to the 3.6-|a.m 
measurement (Figure 9, top panel). However, these models 
still cannot explain the observed brightness-temperature dis- 
crepancy with wavelength. 

The study of eccentric hot- Jupiter atmospheric circulation 
by Kataria et al. (2012) hints at a resolution to this discrep- 
ancy. Their Figure 4 (top panel) shows that, as the planet 
passes through periapsis, the time that the peak temperature 
is reached varies as a function of pressure. This is typical of 
their simulations (T. Kataria 2009, private communication). If 
this differential response is significant in WASP-8b, it would 
introduce a discrepancy in the observations since the Spitzer 
bands sample different altitudes (see Figure 9 right panel). 

Another possibility is to compare the radiative and advec- 
tive timescales at the altitudes sampled by each band. Evalu- 
ating equation (1) of Fortney et al. (2008) using WASP-8b's 
P-T profile, indicates that t^..^^ increases with depth between 
0. 1 and 1 .0 bar, so there should be less longitudinal tempera- 
ture contrast at depth. On the other hand, models of Kataria 
et al. (2012) show that wind speeds decrease with depth, and 
thus Tgjj^, also increases with depth. If the increase of r^^j^, 
with depth is sharper than that of t^^j, then one would ex- 
pect less-homogenized temperatures at depth (but still above 
the photosphere). Hence, the rise in temperature (due to the 
increasing incident irradiation) near periapsis could be more 
pronounced at 3.6 ixm than at longer wavelengths, given the 
weighting functions of Figure 8. 

7. CONCLUSIONS 

S/j/teer observed secondary eclipses of WASP-8b in the 3.6, 
4.5, and 8.0 [a.m IRAC wavebands. In our joint-fit model, 
we estimate eclipse depths of 0.1 13% ± 0.018%, 0.069% ± 
0.007%, and 0.093% ± 0.023% at 3.6, 4.5, and 8.0 i^m, re- 
spectively. These depths correspond to brightness tempera- 
tures of 1552, 1131, and 938 K, respectively. Although the 
3.6-|^m eclipse depth is unexpectedly large, most of the ramp 
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models had consistent depths (within la), while those with 
inconsistent depths fit the data poorly. 

Considering the P-T profile of WASP-8b, KCl, ZnS, Lij, 
LiF, or Na2S clouds could form (see Figure 2(a) of Lodders & 
Fegley 2006). In analogy to brown dwarfs, partial cloud cov- 
erage can cause photometric variability (Artigau et al. 2009); 
however, our three 4.5 i.im observations, spanning 1.5 years, 
have consistent eclipse depths, suggesting no temporal varia- 
tion at secondary eclipse above a hemispheric-mean level of 
^ 35 K (Icr). A moderate cloud layer at altitudes higher than 
those probed by Spitzer would produce a featureless planetary 
spectrum at wavelengths shorter than 2 |a.m (Pont et al. 2008; 
Miller-Ricci Kempton et al. 2012) and would block some of 
the stellar flux, decreasing the temperatures at levels probed 
by Spitzer. Yet, the observed temperatures, which exceed the 
time-averaged equilibrium temperature, challenge this idea. 

Given the high eccentricity, spin-orbit misalignment, and 
observed RV drift in the of WASP-8 system, Queloz et al. 
(2010) suggested the existence of an additional, unseen body 
in the system. Our orbital analysis is consistent with theirs. It 
also improves the orbital parameters and extends the baseline 
of sampled epochs. This constrains the long-term evolution of 
the orbit and aids the search for a second planet, for example 
through the study of timing variations (e.g., Agol et al. 2005). 

The eclipse depths probe the day side atmosphere of WASP- 
8b. Our results rule out the presence of a thermal inversion 
layer, as expected, given the irradiation level from the host 
star. A model with solar-abundance composition explains the 
4.5 and 8.0 [im planet-star flux contrast; however, includ- 
ing the high 3.6 [xm flux contrast requires models that out- 
put nearly 20% of the orbit-averaged incident irradiation, in- 
dependent of the atmospheric composition. If the orbit were 
circular (and thus the irradiation steady-state), the high bright- 
ness temperatures would indicate a very low energy redistri- 
bution to the night side of the planet. For an eccentric planet, 
it at least indicates a short t^^j (Figure 9). 

By modeling the orbital thermal variations due to the eccen- 
tricity of the orbit, we determined that it is possible for WASP- 
8b to achieve temperatures as high as the 3.6 [a.m bright- 
ness temperature. However, the differing brightness temper- 
atures in the other two bands remain puzzling. Neither the 
radiative-transfer model (Section 5) nor the phase-variation 
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model (Section 6) embraces all the physics of the problem. 
The radiative transfer code is a ID, steady-state model repre- 
senting typical dayside conditions. The phase-variation model 
describes emission as a blackbody on a single-layer grid; it 
does not consider absorption or emission features from the 
species in the atmosphere. Clouds (e.g., Cushing et al. 2008), 
atmospheric dynamics (e.g.. Showman et al. 2009), and pho- 
tochemistry (e.g., Moses et al. 2011) are not directly consid- 
ered by these models. 

What we can say for certain is that the assumptions of our 
simple models have been violated, which is not surprising for 
this eccentric planet. While it may be possible to construct 
consistent, realistic models, model uniqueness may be elusive 
until more and better data are available. 

Relatively few exoplanets with equilibrium temperatures 
below 1500 K have been observed at secondary eclipse 
(Cowan & Agol 2011b). The same is true for eccentric 
planets. The characterization of WASP-8b in this work thus 
addresses a particularly interesting, and challenging, region 
of the exoplanet phase space. Observation of other planets 
with similar equilibrium temperatures or eccentricities will 
help discover the physics that drive these unusual atmo- 
spheres. 
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